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ABSTRACT 


The intensity wave-like patterns observed in Synthetic Aperture Radar 
(SAR) are known to be caused by two mechanisms: the microwave radar cross- 
sectional amplitude modalation due to tilt and hydrodynamic Interaction of the 
long ocean waves, and intensity modulation due to the motion of the long ocean 
waves. Two-dimensional closed form expressions of intensity wave patterns 
based on ocean wave swell are developed. They Illustrate the relative impor- 
tance of the amplitude and motion modulations; furthermore, they show that 
velocity bunching and a distortion due to the phase velocity of the ocean wave 
field are independent of the focus adjustment, provided that the second-order 
temporal effects are neglected. Second-order effects are small only over a 
limited range of ocean/radar parameters. Future modeling work should con- 
centrate on two-dimensional expressions and numerical methods, including the 
anisotropy of the amplitude modulation, that will allow quantitatively com- 
pared with measurements. 


1 . INTRODUCTION 


The Interpretation of ocean wave imagery is an area of active research In 
microwave remote sensing of the ocean surface. The principles of Synthetic 
Aperture Radar (SAR) are well understood for point targets; see, for example, 
Raney (1971). However, SAR signatures caused by random and systematic motion 
of an extended surface such as the ocean are not as well understood but are of 
considerable interest. 

To construct a SAR image, the radar utilizes the Doppler shift or equiva- 
lently its phase history, produced by the uniform platform velocity, to locate 
targets in the flight direction. If the targets move during the time interval 
required to form the phase history, then the history is modified and target 
locations differ from the ones expected for stationary targets. Suppose that 
an ensemble of targets are moving uniformly in the flight direction; then 
there is no relative position error between the targets. However, suppose 
that targets exist that are spatially separated in the flight direction wJth 
different radial velocity components. The radar then senses the displaced 
Doppler histories in the flight direction, by a process called “velocity 
bunching." This bunching, which is unique to the SAR, can allow the detection 
of ocean waves and ocean-current boundaries even when the radar cross section 
is uniform. However, second-order temporal (quadratic phase) effects that 
Include the random nature of the surface also called scene coherence and the 
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orbital acceleration of the long waves can degrade the wave imaging process 
caused by velocity bunching. 

The previous interpretation of ocean wave imagery has usually emphasized 
the image formation process along the flight direction, while others have con- 
sidered the two orthogonal directions, along the flight and cross-flight slant 
range directions (see, e.g., Alpers and Rufenach (1979); Jala (1981); Valen- 
zuela (1980); Harger ( 1980); Alpers et al. (198*7. The pin pose of the present 
work is to extend earlier results to ocean waves traveling n an arbitrary 
direction and include the motion of temporal amplitude modulation. To 
accomplish this generalization, it is necessary to develop a two-dimensional, 
analytical expression based on the two mechanisms responsible for wave-like 
patterns in the imagery: (1) the cross-sectional modulation due to tilt and 
.lydrodynamic effects, also called amplitude modulation or the modulation 
transfer function, and (2) the intensity modulation du* to ocean wave motions. 
The relative importance of amplitude and artificial modulation based on an 
arbitrary long ocean wave orientation must be included for quantitative 
modeling. Furthermore, for typical ocean/radar parameters, the image for- 
mation process is non-linear except for limited ocean-wave parameters. 
Therefore, a closed form expression for this mapping is not usually available 
except for the case when the mapping is linear. However, a two-dimensional 
description is useful because it gives physical insight even under the above 
restrictions. 

Jain (1981), and Shuchman and Zelenka (1978) have made specialized mea- 
surements where the maximum image contrast occurs at a focus adjustment dif- 
ferent than expected for a stationary surface. SAR images are brought into 
focus at the focal (reference) plane by adjusting the matched filter chirp 
rates such that Ab » Ab f * 0 in Eq. (6). Jain (1981) has suggested that this 
defocus is equal to the azimuthal component of the long wave phase velocity; 
whereas Alpers et al. (1981) has suggested that it is the radial component of 
the orbital acceleration of long ocean waves that causes this defocus. The 
expressions developed in the present work may help resolve the focus adjust- 
ment discrepancy. Also, we show that the amplitude modulation, which is not a 
focusing phenomenon, can also cause distortions in the inferred ocean wave 
field which may be important for slow flying aircraft measurements. Further- 
more, the two-dimensional encountered wave-like patterns are given in terms of 
the apparent ocean wavelength and direction caused by the phase velocity of 
ocean wave swell which is relevant for both synthetic and real aperture 
radars. 


2. SYSTEM MODEL 

Suppose that a radar platform is moving with velocity V along the x direc- 
tion. Furthermore, assume that a point radar scatterer is located at a range 
r on the surface at x, y, and at a range R at midbeam, x » Vt, as illustrated 
in Fig. 1. The surface elevation associated with the wave field is 


170 



43 



Figure 1. Artist's concept of radar/ocean geometry; 0 is the angle of inci- 
dence, ♦ is the azimuth angle between k the long oceai. vee'" ^ve number 
and the flight direction and 0 is antenna beamwldth. a it v->f .'he 

long waves is traveling in che same direction an the f . rj : erection; 

i.e., « * 90*. 


c(*.y.t) - I r co.(k Kil x * k y «„«> (1) 

n-1 J 

where the ocean wavenumbers k xn> k are the components along the x- and 

respectively, k - k cos ♦ , k - k sin ♦ , and k ■ 

xn n n y n n n 

The complex amplitude signal bsckacattered from the ocean surface depends 
on the reflectivity properties of the surfsce including the time dependence of 


t 

•; 
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the amplitude modulation due to the long ocean waves. The complex amplitude 
received at the platform based on a point scatterer, following Hasselmann 
(1980), is 

-i(b t 2 + b t 2 ) 

dA(t , t) - G( x - Vt,y) E(t-t’) e r r a(x,y,t, dxdy (2) 

where G is a factor that includes the antenna-weighting function E(t) is the 
pulse-weighting function, b is the azimuth chirp rate, b is the range chirp 
rate, and a(x,y,t) is the complex amplitude reflectivity. The dependence of 
A(t,t) on t (fast time) is used in processing to determine the slant range-r 
positions of the scattering element; the Doppler signature is related to the 
platform velocity and slow variability of the reflective surface through t 
(slow time) to infer the aximuthal x-posltlon of the scattering element. 

The received signal is competed with a reference signal in the processor 
which is also called a matched filter repreaenteu here by h(t,T>. The two- 
dimensional convolution whose output is the image complex amplitude for a 
point scatterer is 

A( t-t ' , t-t') - // A( t j-t t^-t ) h(t'-t 1 , t'-Tj) dtjdtj • (3) 

Ba.xacattered microwave signals from the rough sea surface are described 
by a two-scale Bragg scattering model first Introduced by Wright (1968) and 
Bass et al. (1968). This model is consistent with wavelike amplitude patterns 
observed in imaging radar, provided that the radar resolution is smaller than 
the long ocean wavelength. Initially, this amplitude modulation was explained 
by the geometric tilting of the long ocean waves in a local reference plane. 
Some years later, after more complete analysis, it was shown that the strain- 
ing of the short (say 1-100 cm. waves by the long waves &^so influenced the 
modulation. This straining, also called hydrodynamic interaction, causes an 
asymmetrical distribution of the short waves with respect to the long waves 
(Keller and Wright, 1975). 

An extension of the previous S>R results, which shows that the amplitude 
modulation is dependent on the match filtering, la more easily accomplished by 
restricting the model to ocean wave swell, which to a first approximation can 
be represented by a sinusoidal long ocean wave. This simplifies Eq. (1) to 
one Fourier component, c(x,y,t) - c o coedt^x ► fc^y - £t). Therefore the 

complex amplitude reflectivity can he approximated by 

a(x,y,t) a [l ♦ y coe(k x x ♦ k y y - «t)] e’ 1 ^*****) O) 


where a is the complex amplitude modulation Index caused by the long ocean 
waves, $ la the phase perturbations caused by the orbltsl motion, and the hat 
" • “ indicates long ocean wave parameter. For the present model the scene 
coherence Is neglected. The reader Interested In this coherence Is referred 
to Raney (1980), Rufenach end Alpers (1981), and Lysengs and Shuchman (1983). 
The phase perturbation can usually be approximated by 
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where k is the radar wavenumber, and u^ and a^ are - dial rnmponent« r> f the 
orbital ;lociiiy and accel ation of the long waves, r espectivel y. The ampli- 
tude modulation index can be approximated by m =» m ^ 2 + C based on ocean 

# o o x y 

wave swell where m is the modulation transfer function given by others (see, 
o 

e.g., Keller and Wright, 1975); ana dependent on ocean/radar parameters. For 
example, suppose the long ocean waves are traveling along the radar look 
direction (range traveling waves) and the receiving and transmitting antennas 
are horizontally polarized then m^ * 4 (tantf + et: 9) or vertically polarized 

m * -2 3 in 29/(1 + sin 1 2 0) + 4 ctnO (Rufer.ach et al., 1983). 
o 


The image is brought into focus in a reference plane by adjusting the 
quadratic phase of the matched filter. In practice, this focus adjustment is 
usually achieved on an optical prc isor by adjustments of the focal length 
between the platform and a stationary reference surface. Therefore the 

matched filter is given by 


* 2 *2 

-i(b t + b t ) 
h( t , t) - e x r 


( 6 ) 


where .he referer.ee chirp rates b ' and b ' nay differ from the received i.ci r ,> 

rates b ■ k V 2 /R and b ■ 2B /'I by rhe differential rate Ab - ’> ’ - b and 
xo rrr xxx 

Ab f ■ b f ' - b r . The received envelope of the chirp signal is related to the 

azimuthal bandwidth B - b T /2 and range bandwidth B * b T /2 wher n T and 

xxx rrr x 

T f are the azimuth integration time and the rf pulse duration, respectively. 


The intensity modulation based on a point scatter.r riding an ocean swell 
and assuming a Gaussian antenna weighting 


G(t , t) 




(7) 


1. obtained by substituting Eqs. {2), (4), and i 6 ) into Eq. (3) based on the 
high radar wavenumber limit (k ♦ «•), and assuming m « 1. Following Rurenac'.i 
and Alpers (1981), and Hassefmann (1980) for temprral-to-epatial coordinate 

transformation and using straightforward but tedious mathematics. 
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;a | 2 ■ f (1 + £s)(l + T 2 T 2 [1 * i co.(k’« - k' r r)] 

x r 


««p[- (f-) 2 (1 * v-) 2 (*’- * 7 »;i 2 - (^-) 2 (l ♦ ^Hr’-r) 2 ] (8) 


where 


k x =■ k x (l - — ) is the apparent long ocean wavenumber. 


and 


c io : j velocity of light 
jn V_ 

P x * 2 B ’ 


P r - » ~ 


B * 
r 


and 


*> ? 
u 


1 + 


Ab * 
x 


The high wavenumber limit is taken with the azimuth antenna diameter constant 
which means that the integration time is sufficiently small that all quadratic 
phase effects such as orbital acceleration can be neglected. The above 

expressions apply for an image whose azimuth scale is equal to the radar 
antenna beamwidth. Equation (8) illustrates that the amplitude modulation 
field as measured in radar images can be distorted from the actual wave field 
on the ocean surface. This distortion occurs only along the flight direction. 
Furthermore, this distortion is not caused by a quadratic phase effect. In- 
deed, the type of distortion in Eq. (8) is most important for slow flying 
aircraft when the long ocean waves are traveling along the flight direction, 
relevant for both synthetic aperture and real aperture radars. Equation (8) 
shows the dependence of amplitude modulation on image formation based on a 
point scatterer riding on a monochromatic ocean wave swell. 

The image response based on a distributed surface such as the ocean is 
obtained by summing all the elementary point scatterer contributions: 

I(x' t r f ) - // | A | 2 dxdr . (9) 


Using the assumptions of Alpers and Rufenach (1979) namely 
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* R ^ U r 

ra, « r and |l + -~| , 


not varying much within a resolution cell we obtain 
n/4 T 2 T 2 Ab Ab 

I(*',r') ■ (l + + t— ^-) fl + m cos (k x ' + k r ' ) 1 (10) 

. v b J K b J l sar v x r J J 

du X r 

l‘ + ?srl 


• * r 1 /■ 

where m * m exp k 

sar r L , 9 v 

4 n 


l - 2 (k’ 2 


(i +— ' 


Ab , 
x ^ 2 

b ' 
x 


A 9 K r 

k i— vni- 

(' +TT 1 ) 2 


r. ^tl 

R r 

Furthermore, if — ^ — « 1 then Eq. (10) simplifies to 


du 

I(x' ,r') - K(l - y d — )[l + m gar cos(k x x’ + l^r’ )] 


(ID 


Ab 


Ab 


where K = tt/4 T T fl + t— “] fl + r— •) or in terms of ocean wave swell 
x r v D y 

„ x r 

parameters 

lfx',r') - Kfl + n cos fk x' + k r* ) + |m I cos fk x' + k r' + 6)1 (12) 

v > l 0 v x r ^ 1 sar ' ^ x r ; J 


^ ^ / 2 2 2 

where n * R/V ? u> k /sin 0 sin ♦ + cos 0 is the amplitude of the velocity 

o OX i | 

i • i A 2 * 2 

bunching, m H m r /k + k provided the radar resolution filtering of 
1 sar o o x y 

the long waves is neglected, and 6 is the phase angle of the amplitude modula- 
tion. For tilt modulation 5 - »/2, which means that the surface elevation is 
90° out of phase with the amplitude modulation. Equations (10) and (11) show 
that wave-like patterns in the image caused by ocean wave swell are due to two 
mechanisms: velocity bunching and amplitude modulation. These equations hold 
only ove.' a narrow range of ocean/radar parameters since quadratic effects 
such as orultal acceleration have been neglected. However, they do illustrate 
the relative importance of the real and artificial modulations. The amplitude 
modulation is the factor within the square brackets in Eqs. (10) and (11). 
Equation (11) and (12) are linear expressions relating the wave field to the 
image intensity modulation. 

The Synthetic Aperture Radar is generally considered to be a sensor that 
can measure the dominant ocean wavelength and direction. Under certain 
limited renditions, significant distortions can occur in both the observed 
wavelengt.i and direction. Using Eq. (10) or (11) gives 
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where primed parameters indicate apparent parameters , whereas the unprimed 
parameters refer to the actual wavefield. Therefore, the apparent wavelength 
is 



where 


X 


2 ir/kL 1 and k 


k 
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cos$ * 


and the apparent direction is 


(14) 


-1 f sin$ 
tan [ 


-=] 


cos4> + 


to/k 


05) 


A 

As an example consider two cases in which deep water swell with X * 250 m is 
traveling along and opposite to the sensor flight direction. Suppose further 
that the sejisor is a slow-flying aircraft with V ■ 100 m/s. The phase velo- 
city is co/k x = 20 m/s which gives 2 250 (1 ± 0.2) ■ 200, 300 m. This 

illustrates that significant distortion can be due to other than non-linear 
effects under certain conditions. However, it should be noted that non- 
linearities, caused by orbital acceleration are likely to dominate on most 
occasions. 


3. DISCUSSION AND SUMMARY 

The wave-like patterns observed in SAR imagery are caused by two mecha- 
nisms: (1) the radar cross-sectional (amplitude) modulation due to tilt and 
hydrodynamic modulation by the long ocean waves and (2) intensity modulation 
due to the motion of the ocean surface which is unique to SAR. The motion- 
induced modulation can be separated into modulation enhancement and degrada- 
tion due to the systematic orbital acceleration of the long waves and the 
degradation of the modulation due to the stochastic character of the wind 
waves. The radar amplitude modulation is dominant for ocean waves traveling 
perpendicular to the flight direction whereas the motion-induced modulation 
may dominate when ocean waves are traveling along the flight direction. 

Determination of the relative importance of these underlying mechanisms is 
essential to a complete understanding of the radar interaction with the ocean 
waves and the modeling of SAR wavelike patterns; for example, see, Lyzenga et 
al., 1984. Interest in the different mechanisms has increased in recent years 
since it is difficult to separate them except for the special case when the 
waves are traveling exactly perpendicular to the flight direction (range 

176 



P.V 'l 



waves). For this case, motions on the surface are negligible. Most of the 
theoretical work has emphasized the modulation due to ocean wave motion and 
its associated non-linearities. A two-dimensional closed form expression is 
not available except over a limited range of ocean/radar parameters. This has 
led some workers to consider Monte Carlo numerical methods to model the image- 
formation process (Alpers, 1983). Two-dimensional expressions have been deve- 
loped which are valid over a limited range of ocean/radar parameters. Indeed 
in this limited range, the Monte Carlo results that required numerical com- 
putation could be compared with closed form results such as Eqs. (10) and 
( 11 ). 

Equations (10) and (11) Include the effects of both radar amplitude end 
velocity bunching modulation. These equations show the following salient 
features: (1) Velocity bunching is Independent of focus adjustment, provided 
that the quadratic phase is negligibly small which normally means short 
integration times; however, focus adjustment will degrade the image in the 
same manner as in a stationary scene, and (2) Ocean wave field distortion 
caused by the motion of the wave field (phase velocity of the long waves) 
relative to the platform velocity is Independent of focus adjustment, again 
provided that the quadratic phase is negligibly small. These conclusions are 
in disagreement with Jain (1981), who claims that the focus adjustment is 
equal to the azimuth component of the long wave phase velocity. Furthermore, 
the dependence of the focus adjustment on orbital acceleration as suggested by 
Alpers et al. (1981), is not relevant since these expressions are not valid 
when the orbital acceleration is Important. 

It is recommended that future work continue to emphasize development of 
quantitative two-dimensional models which give the relative Importance of 
amplitude and motion modulation for any arbitrary orientation of the ocean 
wave field. Furthermore, the anisotropy of the amplitude modulation must be 
included in order to quantify this modulation. 
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